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The Gaussian point transmission model for calculating optical 
fiber splice loss is extended to the general case of splice loss between 
fibers which differ in one or more intrinsic parameters — core radius, 
index of refraction profile shape, and maximum index of refraction 
difference between core and cladding. The model is first verified for 
splices with index- of -refraction profile mismatch. The average differ- 
ence between calculated and measured splice loss due to profile 
parameter mismatch is 0.04 dB. Comparisons are also made between 
calculated and measured splice loss for ten different splices with 
mismatch in all three intrinsic parameters. The average difference 
between the calculated loss and the average of several measured 
losses for these ten cases was 0.06 dB. The additional losses intro- 
duced by transverse offset measured for one set of mismatched fiber 
splices agree with the calculated values within 0.1 dB. Loss due to 
misalignment of elliptical core fibers is calculated and measured 
with agreement within 0.06 dB for the maximum loss case. Both the 
model and the experimental data show that, for a given percentage 
mismatch, index of refraction profile parameter mismatch and core 
ellipticity contribute significantly less to splice loss than mismatch of 
core radius or numerical aperture. A family of curves for splice 
loss vs transverse offset is presented for various numerical aperture 
mismatches and core radius mismatches, since these parameters are 
typically the largest components of splice loss in practical fiber optic 
systems. 

I. INTRODUCTION 

One factor which must be considered in the development of fiber 
optical communication systems is the effect of fiber core parameter 
manufacturing variations on splice loss. These intrinsic fiber core 
parameters 1 are the maximum index-of-refraction difference between 
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core and cladding, A, the index of refraction profile parameter, a, the 
radius, R, and core ellipticity, e. 

The development of a phenomenological Gaussian point loss model 2 
allows an approximate analytic treatment of the loss induced by a 
butt-joint splice between fibers which differ in one or more intrinsic 
parameters. Previous models based on the uniform power distribution 
assumption 1 have exhibited only limited agreement between calculated 
and measured splice loss, 3,4 whereas the Gaussian model gives good 
agreement with experimental data. 

Previous work 2 developed the Gaussian point loss model and gave 
theoretical and experimental results for A mismatch, R mismatch and 
transverse offset. This model has been used to estimate system losses 
due to random splicing between fibers whose intrinsic parameter 
variation distributions are known. 5 An extension of this model to the 
case of index-of-refraction profile (a) mismatch or ellipticity (e), and 
to combinations of intrinsic factors plus transverse offset, is presented 
in this paper. Although the extension of the model is straightforward, 
the resulting calculations require the use of approximate numerical 
quadrature techniques applied over irregular areas of integration. This 
paper presents a unified analysis of the various effects of intrinsic 
factors and transverse offset on splice loss utilizing the Gaussian model. 
Representative experimental data are also presented. 

II. THEORETICAL DEVELOPMENT 

An approximate analytical treatment of the loss in a butt-joint splice 
between two fibers with differing intrinsic parameters can be carried 
out using a Gaussian model. 2 This model assumes that a steady-state 
power distribution (after a long length of fiber) can be modeled as a 
Gaussian distribution of the power within the solid angle defined by 
the local numerical aperture, na, at any point on the fiber core 
(Fig. 1). With this assumed power distribution at the input to a splice, 
the ratio of the power received, P2, at any radial position, r, to the 
power transmitted, p\, at that point is related to the ratio of the squares 
of the na's at that point by the following equation: 2 

,l+-r^4po-e (NA2)2/(NA,)2,np " for na 2 < na, 

Pl-\ (NA,) 2 ^ (1) 




for NA 2 > NAi, 

where po is the point defining the width of the Gaussian distribution 
corresponding to 0.1 of its maximum. The expression t(r) is called the 
point transmission function. 
Using the usual class of circularly symmetric profiles, 1 
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p(r)=EXP 



Fig. 1 — Gaussian power distribution. 



NA(r) as n ^2A 



'-i 



1/2 



where A = (n - /i c )/no is small 

no = maximum core refractive index 
n c = refractive index of cladding 
a = index profile parameter [a(r) = a] 
R = fiber core radius. 
Substituting eq. (2) into eq. (1) and defining 



Q = 



(NA 2 ) 5 



A 2 (l - r" 2 ) 



(2) 



(3) 



(na,) 2 A,(l - ^"r"') 

gives the transmission coefficient at a radial distance r from the core 
axis, 

1 + Qpo - p? for Q < 1 
1. for Q > 1, 



<(r) = 



(4) 



where k = Ri/R\ with /? 2 normalized to 1. Subscripts 1 and 2 refer to 
the transmitting and receiving fibers respectively. 

To find the total transmission through the splice, the point trans- 
mission function is multiplied by the input power of the transmitting 
fiber at each point and integrated over the area of core overlap. A 
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transmission ratio is obtained by dividing by the total input power: 

-2w r x 



Pi 



JO 



t(r)(l - k">r"*) 2 rdrdd 



r = y r ^jrr k , (5) 

1 ' (1 - k a *r a *) 2 rdrd6 



JO 



where x = the lesser of 1 or 1/k. 

For circularly symmetric profiles with no transverse offset: 



t(r)(l - k a 'r"Trdr 

Px ai/2k 2 (a, + l)(a, +2)' 

This formula represents the total intrinsic loss immediately after the 
splice. Substitution of eq. (4) for t(r) into (6) results in an integral 
which must be solved numerically for t{r) ¥= 1 except for the special 
case of a\ = a* and R\ = Rz- 

The inclusion of transverse offset in the problem requires the eval- 
uation of the double integral in the numerator of eq. (5) over the area 
of overlap of the two fiber cores as shown in Fig. 2: 

t(r)(l - k a 'r a >) 2 rdrdd 

(7) 



P, 27rai/2k 2 (a i + l)(a, + 2)* 

Equation (7), when integrated over the area of core overlap, repre- 
sents a general solution for any combination of intrinsic mismatch and 
transverse offset for short receiving fibers. Long receiving fibers require 
the use of the long length t(r) as given in Ref. 2. The splice loss can be 
calculated by evaluating Q at each point and integrating eq. (7) 
numerically using the appropriate expression for t(r). This numerical 
integration can sometimes be simplified by careful consideration of the 
particular splice parameters. The absence of transverse offset, for 
example, reduces the numerator to a single integral in r for circularly 
symmetric fibers. These results assume flat, clean, perpendicular fiber 
ends with index-of-refraction matching fluid and no angular or longi- 
tudinal offset. In addition, a steady-state distribution is assumed at 
the input to the splice, requiring either a long input fiber or a launched 
power distribution approximating the steady state. 

III. COMPARISON OF EXPERIMENTAL AND THEORETICAL VALUES 

Experimental verification of the Gaussian model has followed a step- 
by-step procedure. Splice losses were measured between selected sets 
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Fig. 2 — Overlap regions for offset fiber cores. 

of fibers which differed primarily in one intrinsic parameter. The 
results for A and radius mismatch were previously reported in Ref. 2, 
along with transverse offset results for both short and long receiving 
fibers. Results for a mismatch and combinations of intrinsic mismatch 
and transverse offset are given here. 

All measurements were made using a pulsed, 0.82-jum, laser source 
with a pigtail which was loose-tube-spliced 6 to the transmitting fiber. 
The minimum length of the transmitting fiber was 550 m, except for 
one step index fiber (~250 m), to provide an approximate steady-state 
power distribution at the splice. The transmitting fiber was wound 
under tension to reproduce typical microbending losses found in some 
cables to simulate transmitting power distributions of interest. The 
receiving fiber was approximately 1 meter in length, though two tests 
were repeated using 10-m lengths to determine if any cladding modes 
were present. The results were essentially the same in either case. 

IV. a MISMATCH 

Splice loss due to a mismatch was measured for a few representative 
fiber pairs. The results are shown in Fig. 3 along with the general 
a mismatch curves generated by this model. The agreement between 
theory and measurement is good. For a-i = 1.5 and a\ = 2.0, a 25- 
percent mismatch in a, the splice loss is less than 0.2 dB. Sensitivity of 
splice loss to a mismatch is therefore substantially less than that for 
A or core radius mismatch. 2 

Fibers with combinations of a and A or a and radius mismatch were 
also spliced to observe the effects of combinations of factors. The 
results are given in Table I, where the first three columns indicate the 
percentage difference of each intrinsic parameter of the transmitting 
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Fig. 3 — Transmission vs a mismatch. 

Table I — Calculated and measured splice losses between 
selected fiber pairs for a mismatch verification 





Percentage 
Difference* 




Predicted 
Loss 

(dB) 


Measured Loss (dB) 


Number 


a 


A 


R 


Low 


Avg. 


High 


Readings 


+26.1 
+20.9 
-32.5 
+ 13.1 

+ 00 




-3.5 

+37.6 


-3.5 


+6.3 




-1.1 

-14.8 


0..374 
0.105 
0.558 
0.074 
0.864 


0.258 
0.053 
0.564 
0.101 
0.442 


0.357 
0.110 
0.614 
0.189 
0.585 


0.458 
0.156 
0.667 
0.276 
0.674 


6 
7 
5 
7 
8 



* +26.1 indicates that the a of the transmitting fiber is 26.1 percent higher than the 
a of the receiving fiber. 



fiber from those of the receiving fiber with the sign indicating the 
direction of the change (+ indicates a larger transmitting fiber param- 
eter). The average measured loss is given, as well as the range and 
number of measurements. Each measurement involves disassembly of 
the splice, fracture of new ends on both fibers, and reassembly of the 
splice. 

The range of loss measurements for some combinations of parameter 
mismatch is relatively large. Several factors contributed to this varia- 
tion, (i) The detector used for early measurements degraded in its 
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repeatability properties and was replaced by a different detector for 
later measurements. These repeatability problems were primarily due 
to lack of precision in positioning the fiber relative to the detector, 
(w) Failure of the lock-in amplifier to track frequency drift in the 
source was found to cause measurement variations of as much as 
±0.1 dB. (Hi) End quality and contamination, along with small varia- 
tions in the positions of the fibers in the splices, also contribute to 
measurement variations, (iv) Fiber profile asymmetries may also con- 
tribute slightly to measurement variations. 

An attempt was made to test the limits of applicability of the model 
by applying it to a splice between step-index (a = <») and approximately 
parabolic-index fibers. The results, as shown in Table I, indicate that 
the model gives a fair estimate of the loss even in this extreme case 
when the A and radius mismatches are included in the calculation. The 
shortness of the step-index transmitting fiber in this case (250 m) may 
have caused an incomplete filling of the transmitting mode structure. 

The results for general intrinsic parameter mismatch measurements 
are given in Table II. It can be seen that some measurements agree 
very well with the calculated values, while others vary on both the 
high and low sides. The average relative difference bet ween theoretic al 
and experimental values, (\/n) £ (calculated loss - measured loss), 
was 0.03 dB for all graded-index fiber data in Tables I and II. 
The average absolute difference, (1/n) £ | calculated loss - 
measured loss|, was 0.11 dB. It should be noted that inaccuracies in 
measuring a, A, and core radius as well as variations of a over the core 
probably account for some of the difference between theoretical and 
experimental values given here. The experimental data were accumu- 
lated over a three-month period using several different detectors. All 
the data are reported here. No systematic errors are apparent. 

V. CORE ELLIPTICITY 

Measurements of splice loss due to core ellipticity were made as a 
function of the angle between the major axes of the transmitting and 
receiving fibers (Fig. 4). This produced an approximately sinusoidal 
curve with zero loss at the 0° and 180° points and maximum loss at 
the 90° and 270° points. This maximum can be compared to the value 
calculated from the model. 

The measurements were limited by a lack of suitable fibers with 
significant core ellipticity as well as by the low magnitude of the losses 
involved. Several fibers were selected which had core ellipticities of 4, 
10, and 16.5 percent. The fibers with 4- and 10-percent ellipticity were 
too short to permit measurements with long fiber lengths after the 
splice. The maximum loss measured was less than 0.05 dB for the 
4-percent fiber and less than 0.1 dB for the 10-percent fiber. The 
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Table II — Calculated and measured splice loss for 
general parameter mismatch verification 





Percentage 
Difference 




Predicted 


Measured Loss (dB) 


Number 








Loss 
























a 


A 


R 


(dB) 


Low 


Avg. 


High 


Readings 


-5.1 


-0.22 


+0.4 


0.001 


0.059 


0.090 


0.156 


6 


+5.1 


+0..22 


-0.4 


0.032 


0.102 


0.168 


0.216 


5 


+33.3 


-1.2 


+7.8 


0.495 


0.468 


0.497 


0.537 


5 


+4.5 


+7.6 


+2.7 


0.195 


0.075 


0.155 


0.321 


8 


+24.6 


+ 19.9 


+0.8 


0.579 


0.275 


0.296 


0.310 


5 


-10.5 


+50.0 


+0.8 


1.27 


0.817 


1.03 


1.42 


24 


+21.0 


+ 13.3 


-1.9 


0.340 


0.171 


0.216 


0.260 


5 


+31.3 


+ 13.3 


-3.0 


0.437 


0.237 


0.264 


0.321 


5 


-14.6 


+45.9 


-1.9 


0.996 


0.703 


0.876 


1.06 


12 


-30.4 


+36.8 


+0.8 


0.536 


0.610 


0.684 


0.774 


7 



RECEIVING FIBER 




-TRANSMITTING FIBER 



(a) 




(b) 

Fig. 4— Elliptical core fibers, (a) 90-degree core misalignment, (b) 16.5 percent 
ellipticity fiber core. 
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repeatability of the measurements is a few hundredths of a decibel. 
This prevented an accurate determination of the effect of the elliptic- 
ity. These losses were approximately the same as those predicted by 
the Gaussian model (0.025 dB for 4 percent and 0.08 dB for 10 percent) 
(Fig. 5) and significantly less than those predicted by the uniform 
power model (0.12 dB for 4 percent and 0.325 dB for 10 percent). 
Figure 4b shows the core cross section of the 16.5-percent fiber. 
Although not a perfect ellipse, this was the only fiber available with 
sufficient ellipticity to give a splice loss large enough to be accurately 
measured. 

The 16.5-percent ellipticity fiber was an extremely high loss fiber, 
~50 dB/km for long sections of fiber. The near-field patterns for a 
100-m length and an 800-m length were essentially the same, and it 
appeared that the high-loss characteristics of the fiber produced full 
mode coupling in the 100-m length. For this reason, it was felt that the 




,a-b 



12 16 

P" : '■ — x 100) IN PERCENT 

a 

Fig. 5 — Maximum splice loss due to core ellipticity. 
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100-m length was sufficient to satisfy the assumption (in both Gaussian 
and uniform power models) of a steady-state power distribution enter- 
ing the splice. 

Figure 6 shows the splice loss for this fiber with a 1-m length after 
the splice as a function of the angle between the semi-major axes of 
the ellipses. The general shape of a curve through the data points and 
the average maximum loss of 0.13 dB are in good agreement with the 
expected shape and the value calculated from the model of 0.19 dB. 
The discrepancies may be partly attributable to the deviation from 
pure ellipticity of the core cross section as shown in Fig. 4b and to the 
shortness of the transmitting fiber. The difference between the average 
maximum measured loss and the loss calculated with the Gaussian 
model was only 0.06 dB. This was approximately the same as the noise 
and repeatability properties of the measurement system. 

Figure 7 shows the results when a 100-m length of receiving fiber is 
used. Again the general shape of the curve agrees with that expected, 
although more noise is apparent in the data. The average maximum 
splice loss was 0.18 dB which is in fair agreement with the calculated 
value of 0.31 dB. Measurements with longer lengths of receiving fiber 
were attempted, but the results were inconclusive because of the high 
attenuation and reduced signal-to-noise ratio. 



VI. INTRINSIC MISMATCH PLUS TRANSVERSE OFFSET 

Experimental and theoretical values for splice loss vs transverse 
offset with intrinsic parameter mismatch are shown in Fig. 8. The 
fibers used in this part correspond to line 9 in Table II. The experi- 
mental points represent the average of four readings, two left of center 
and two right of center, as the receiving fiber was traversed across the 
transmitting fiber. The agreement is very good. 
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Fig. 6 — Core ellipticity splice loss, short length. 
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Fig. 7— Core ellipticity splice loss, long length. 
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Fig. 8 — Intrinsic parameter mismatch and transverse offset. 

The model can also be used to generate parametric families of curves 
such as Fig. 9 and Fig. 10, showing the importance of the different 
intrinsic parameters when combined with an extrinsic parameter. A 
mismatch and transverse offset are shown in Fig. 9 and radius mis- 
match and transverse offset in Fig. 10, since they are the chief intrinsic 
and extrinsic factors contributing to splice loss at this time. In partic- 
ular, Fig. 9 emphasizes that, if the receiving fiber na is greater than 
the transmitting fiber na, the sensitivity to transverse offset is signifi- 
cantly reduced, especially for small offsets. Sensitivity to transverse 
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Fig. 9 — Splice loss due to A mismatch and transverse offset. 
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Fig. 10 — Splice loss due to radius mismatch and transverse offset. 



0.8 



offset is slightly reduced for a receiving fiber na less than the trans- 
mitting fiber na, in fact, the maximum sensitivity to transverse offset 
occurs for identical fibers. 

VII. CONCLUSION 

It has been demonstrated that the Gaussian point transmission 
model gives good estimates of splice loss due to intrinsic parameter 
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mismatch and transverse offset. The effects of combinations of these 
factors on splice loss has been characterized and shown to be signifi- 
cant. Effects due to a and e mismatch are smaller for a given percentage 
mismatch than for A or core radius mismatches. 
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APPENDIX 

This section provides more detail on the application of this theory 
to various general classes of fiber splices. The integrals in eq. (7) 
cannot be performed exactly, except in a few special cases. Numerical 
evaluation is usually required, but considerable simplification can 
sometimes result from careful consideration of the particular splice 
parameters. The case of intrinsic parameter mismatch with no elliptic- 
ity and no transverse offset is examined first. 

A 1. Splices without transverse offset 

Zero offset splices can be broken down into three classes for the 
particular family of radially symmetric profiles assumed here. 

Class I contains splices with all transmitting parameters greater 
than or less than all corresponding receiving parameters. Class la is 
illustrated in Fig. 11a. 

Class \ A _ . Class 



For la 



For lb the transmission equals unity. 

Class II splices have one intersection or the core profile. Class Ha is 
shown in Fig. lib. 



<, Aj > A 2 T , < Ai < A 2 
Ia: ] R^R 2 Ib - U<* 2 



p 2 = t(r)(l - k°^) 2 rdr. (8) 



Class j Ai > A 2 Class ( A x < A 2 
Ha: \Ri<R 2 lib: \Ri>R 2 
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In both these cases, the transmission is unity where NA2 exceeds NAi 
and a function of r over the remainder of the profile. The crossover 
point, r', is determined by equating the nas and solving the resulting 
nonlinear equation numerically: 

Aid - #V») = A 2 (l - r a *). (9) 

Then 



J r ' 



t(r)(l - k a 'r a ') 2 rdr + (1 - k" l r a ') 2 rdr Ha 

h Jr' 

P 2 =-< (10) 

rVk 

(1 - k ai r° l ) 2 rdr + t(r)(l -k°*r a *) 2 rdr lib 



Class III contains all remaining possibilities. Class Ilia is shown in 
Fig. lie. 

Class ( Al - Az Class ( Al - A2 
< <x\ < a 2 TTTK . \ cti > a 2 . 

IIIa: \r^r 2 mb - Usa 

Class III mismatches may or may not have the crossovers shown in 
Fig. lie. Those that do not cross over correspond to Fig. 11a; however, 
it is necessary to solve eq. (9) to determine whether a crossover occurs. 
If no r' exists, the loss can be calculated using the appropriate formula 
for Class I. If an r' exists, the second crossover point r" must also be 
found to determine the appropriate areas of integration. Again, the 
transmission is unity where NA2 is greater than NAi and is a function 
of r over the rest of the profile. Then 

t(r)(l - k a, r a >) 2 rdr + (1 - k a, r a, ) 2 rdr 

Jo Jr ' Ilia 

+ I t(r)(l - k'"r a> ) 2 rdr 

J r' 



P 2 =- 



I (1 - k ai r a ') 2 rdr + t(r)(l - k a *r ai ) 2 rdr 

JO Jr' 



(ID 



r' 

Illb 

l/A 

+ I (l-tfVOWr 
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Fig. 11— Intrinsic mismatch profiles, (a) Class la. (b) Class Ha. (c) Class Ilia. 

Two special cases also occur when A] = A2 or ct\ = a 2 . 
If A, = A 2 (Class III): 



r ,_M 



(12) 



In this special case, r' exists only if R\ > R2 and «i < a 2 or Ri < R2 and 
oti > cx2. Otherwise, the only crossover is at r = 0. 
If ai = a 2 = a (Class II): 



r = 



A 2 - Ai 



V A 2 — k a ki i 
If «i = a 2 in Class I or III, there is no crossover. 



(13) 



A2. Splices with transverse offset 

The addition of transverse offset greatly complicates the problem. 
The four cases of Ref. 4 must be included and the integrals evaluated 
over the areas of overlap. The power transmitted at any point is the 
product of the point transmission function, t(r), and the power avail- 
able to be transmitted at that point. 



p(r) = t(r)(l - k°^y. 



(14) 



However, since the cores are not centered, one of the na functions 
must be transformed to the center of the other core. The integral can 
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then be performed numerically setting t(r) equal to unity at any point 
at which NA2 is greater than NAi. In the cases in which the two core 
boundaries do not intersect, Cases I and III of Ref. 4, one obtains a 
double integral referenced to the center of the smaller core. This is 
illustrated in Fig. 12a. The definition of d is the normalized separation 
of the fiber core centers. 
Case I: Ri<R 2 ,d<l- 1/k. 

Ti*3.risv6rs6 Offset 

d = - (R 2 = receiving fiber core radius) (15) 

R2 

P 2 = 2 I I fc(r)(l - k a 'r a ') 2 rdrdB (16) 

Jo Jo 



ti(r) = 



1 + QiPo-/??' Qi<l 

(17) 
1 Qi>l 



A 2 [l-(r 2 + rf 2 -2rrfcosflp /2 ] 
Vl AiU-tfV") 

Case III: Ri > R 2 , d < 1/k - 1. 

P 2 = 2 I I t 2 (r)[l - #«'(r 2 + d 2 - 2rdcos O)°* /2 ] 2 rdrd0 (19) 

Jo Jo 

fl + Q 2 po-p$ 2 Q 2 <1 
Mr) = \ . (20) 

U Q 2 >1 



A 2 (l - r" 2 ) 
Ai[l - k a> (r 2 + d 2 - 2rdcos ^) Q,/2 ] 



Q*= >r, ,..,_.2..« :,.^WV (21) 



Cases II and IV of Ref. 4 present two different possible conditions of 
overlap requiring different approaches to the numerical integration. 
These are illustrated in Figs. 12b and 12c. The area of overlap is 
divided into regions which permit the use of the functions t\(r) and 
t 2 (r) defined above. The two separate conditions are distinguished by 
the following criteria: 
Case II: R 2 > Ri and d > 1 - 1/k (interchange R 2 and Ri in Fig. 12). 



1 



Condition 1: 1 - -5 - <r < 0: Fig. 12b 
k 

1 



Condition 2: 1 - -p - d £ > 0: Fig. 12c. 
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(b) (c) 

Fig. 12— Parameter mismatch and transverse offset, (a) Cases I and III offset, 
(b) Condition I 6< > n/2. (c) Condition II 8-. < w/2. 



Case IV: R 2 < Ri and d>l/k- 1. 

1 

k 



Condition 1: -^ - 1 - rf 2 < 0: Fig. 12b 



Condition 2: — 2 - 1 - d 2 > 0: Fig. 12c. 

In practice, the two different conditions in Cases II and IV are 
distinguished by the sign of cos 6 2 . 
Case II: 

1 + d 2 - 1/k 2 



COS 0i = 



COS 02 = 



2d 
1-d 2 - 1/k 2 



2d 



(22) 
(23) 



Condition I: cos 2 < (Fig. 12b). 



'0 J cos 



_ x t 2 (r)[l - k'"(r 2 + d 2 - 2rd cos 8r /2 frdrd0 



l/A 



- I \^_ti(rHl-*r i r at ) 2 rdrdB, 

-0., J cosS 



(24) 
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where 



~^jg^. 



Condition II: cos 2 > (Fig. 12c). 

P 2 = I Lx fe(r)[l - #*»(r 2 + d 2 - 2rd cos 0)"> /2 ] 2 rdrd0 

JO J cos (9 

+ I l*™* ti(r)(l - k°*r a *) 2 rdrd0, 
Jo Jo 

•* r\/k 

+ 1 ti(r)(l - k a ^) 2 rdrdti (26) 



rr 

e 2 Jo 



where 



2d 
Case IV: 



1-d 2 - 1/k 2 

x= 771— ' • ( 27 > 



„ 1/k 2 + d 2 -l 
cos 0i = — — (28) 



2d 

-d 
2d 



c*- 1 '*-*- 1 . (29) 



Condition I: cos 2 < (Fig. 12b). 



P - I L-**i(r)(l - k a *r a >) 2 rdrd8 

'0 •* cosfl 

+ | L_ *(r)[l - #"(r 2 + d 2 - 2rd cos 0) a * /2 frdrdB, (30) 

J it— $2 J cos 

where 
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Condition II: cos B 2 > (Fig. 12c). 
ri/k 

J COB 6 



I ' \ co ° 6 t 2 (r)[l-k a '(r 2 + d 2 -2rd cos 0)°> /2 ] 2 rdrdd 

Jo Jo 

I I fc(r)[l-A a '(r 2 + rf 2 -2rrfcos^r /2 ] 2 rrfrd», (32) 



+ 

le, Jo 



where 



i-V'-f- 1 . <33) 

2d 



Application of the Gaussian model to the case of elliptical mismatch 
is restricted to identical fibers with no transverse offset. The splice 
loss in this case will vary from zero, when the axes of the elliptical 
cores are perfectly aligned, to some maximum, when they are at right 
angles to each other. For identical elliptical fibers, the Gaussian model 
was used to calculate the maximum loss. 

Figure 4a illustrates the geometry of the maximum loss case. For 
identical fibers, ai = a 2 = a, Ai = A 2 and a\ = a 2 and b\ = b 2 where a 
and b are the semimajor and semiminor axes of the elliptical core and 
the subscripts 1 and 2 refer to the transmitting and receiving fibers, 
respectively. The symmetry of the system shown was used to reduce 
the calculation of the fraction of the power transmitted through the 
splice to one quadrant. 

At any angle 0, 



R 2 la 2 sm 2 6 + b 2 cos 2 6 .... 

m -Try *»*»+*** ■ m 

Then Q, the ratio of the numerical apertures squared at any point in 
the splice, is 

«>•.•-££= (35) 

for < < n/4, r < R 2 and for IT/4 < 6 < 11/2, r < R lm The total 
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transmission ratio with R 2 normalized to unity is 

*-8 



/4 /-l p/2 r\/k 

t(r, 0)(1 - k a r a frdrdB + (1 - k a r a ) 2 rdrd0 

JO Jir/4 Jo 



•w/2 r l/k 

(1 - k a r a ) 2 rdrdO 

(36) 



Jo Jo 



The point transmission from 45° to 90° in the first quadrant is unity 
since the receiving numerical aperture is larger than the transmitting 
numerical aperture in this region. Solution of the integrals in eq. (36) 
by numerical methods gives the results shown in Fig. 5 for the maxi- 
mum loss due to core ellipticity. The long length loss was found by 
correcting for the additional loss in a long fiber due to mode redistri- 
bution. 2 

The uniform power model loss was found using eq. (36) with t(r, 0) 
set equal to Q, the ratio of the squares of the numerical apertures, 
using 1 — kfr" as the near-field power function. The difference between 
the Gaussian model and the uniform power model is appreciably 
greater for elliptical mismatch than for the other intrinsic mismatches 
due to the greater difference between the two models near the core- 
cladding boundary. 
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